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In decay modes B^ D^*^^ oi*"^ and B*^ [){*)+£)(.*)- ^ none of quarks in final 
states is tlie same as one of B{Bs) meson. They can occur only via annihilation 
diagrams in the Standard Model. In the heavy quark limit, we try to calculate the 
branching ratios of these decays in perturbative QCD approach without considering 
the soft final state interaction. We found branching ratios of B^ oi*^^ D^*^ are 
at the order of 10~^, and branching ratios of B^ — > _d(*)+Z)(*)^ are of 10^3. Those 
decay modes will be measured in B factories and LHC-b experiments. 

PACS numbers: 13.25.Hw,12.38.Bx 



I. INTRODUCTION 

As an important way in testing the Standard Model and searching for new physics, 
rare B decays become important in particle physics. Although some of them have been 
measured by B factories, many of them are still under study from both experimental and 
theoretical sides. In theoretical side, the factorization approach has been accepted because 
it can explain many decay branching ratios successfully. Recently many efforts have been 
made to explain the reason why the factorization approach has worked well. One of them 
is perturbative QCD approach (PQCD) P|, in which we can calculate the annihilation 
diagrams as well as the factorizable and non-factorizable diagrams. It has been applied 
to exclusive B meson decays, such as 5 — > 7r7r(p) [^, S — > Kn j3|], i? — > Di*^K 0| and 
some other channels jH S 0| • 

Recently, J. O. Eeg et al. computed the B^ DfDj and 5° — > D^D^ decays 
using heavy-light Chiral quark model which is a non-perturbative approach [8] . As shown 
in Fig^ the four quarks in final states and D~ are different from the ones in the B 
meson, and there is no spectator quark. So this decay is a pure annihilation type decay. In 
the factorization approach, this decay is described as h and dm. B meson annihilation into 
vacuum and , D~ being produced from vacuum afterwards. If we calculate this decay 
in factorization approach j9| , we need the D'^ D~ form factor at very large momentum 
transfer 0{m%), but it is zero due to vector current conservation. So it is difficult to 
calculate this decay model in factorization approach. In the so called QCD factorization 
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approach [l^ , the annihilation contribution is plagued by the endpoint singularity. Thus 
it is only parameterized as a free parameter for this kind of contribution. On the other 
hand, by including the transverse momentum of the partons, the PQCD approach is free of 
such singularities. Furthermore, the Sudakov factor induced by the inclusion of transverse 
momentum helps the convergence of factorization. 

PQCD approach has been recently applied to B meson decays with one charmed meson 
in the final states ^, . In the typical B Dn decays, the momentum of the final state 
meson is approximately |mB(l — r^), with r = m|,/m^. This is still large enough to 
make a hard intermediate gluon in the hard part calculation. Therefore the predicted 
results in PQCD agree well with the experimental data. The PQCD calculation of B^ — > 
^{*)+^(*)- _^ ^(*)+^(*)_ ^gpg^yg ^j^j^ charmed mesons in final states may 

be questionable, since the momentum of final state meson here is relatively smaller than 
that of B ^ Dtt case. However, after calculation we find that the momentum of final 
state meson is ^ttibV^ — ~ |mB(l — 2r^), which is only a little smaller than that 
of i? ^ Dn case. For example, the W boson exchange causes bd cc, and the ss quarks 
are produced from a gluon. This gluon attaches to any one of the quarks participating in 
the W boson exchange. In the heavy quark limit, we apply the hierarchy approximation, 
which is adopted by ref. [^, 0|, ^qcd ^ ttid ^ ttib- In this limit, the D meson momentum 
is nearly 7/1^/2. According to the distribution amplitude used in ref.0|, the light quark in 
D meson carrying nearly 40% of the D meson momentum. It is still a collinear quark with 
1 GeV energy, like that in i? ^ Dtt, B ^ tttt decays. Therefore the gluon connecting 
them is a hard gluon, so we can perturbatively treat the process where the four-quark 
operator exchanges a hard gluon with ss quark pair. 

The framework of PQCD and analytic formulas for the decay amplitudes will be shown 
in the next section. In section UTTj we give the numerical results and discussion. Finally, 
we conclude this study in section IIVI 



II. FRAMEWORK 

The factorization theorem allows us to separate the decay amplitude into soft($), 
hard(if), and harder (C) dynamics characterized by different scales It is expressed 

as 

Amplitude- j (fkid'^k2d^h 1i[C{t)^B{ki)^D{k2)^D{h)H{ki,k2,h,t)\, (1) 

where fcj's are momenta of light quarks included in each meson, and Tr is the trace over 
Dirac and color indices. The soft dynamic is factorized into the meson wave function 
$Mj which describes hadronization of the quark and anti-quark pair into the meson M. 
The harder dynamic involves the four quark operators described by the Wilson coefficient 
C{t). It results from the radiative corrections to the four quark operators at short dis- 
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tance. H describes the four quark operator and the quark pair from the sea connected by 
a hard gluon whose scale is at the order of Mb-, so the hard part H can be perturbatively 
calculated. The hard and harder dynamics together make an effective six quark interac- 
tion. The H depends on the specific process, while $m is independent of any processes. 
Therefore we may determine (^m by other well measured channels to make prediction 
here. 






(d) 



FIG. 1: Diagrams for B DfD^ decay. The factorizable diagrams (a) and (b) contribute to 
Fa-, and the nonfactorizable diagrams (c) and (d) do to Ma- 



We consider the B meson at rest for simplicity. It is convenient to use light-cone 
coordinates (p'*', p~, Py), which is defined as: 
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P = ^ . P = ^ , Pt = {p,P)- (2) 

Thus, expanding up to the order of r^, we can take the B meson and two D^*^ meson 
momenta as: 



Mb , , Mb, ■? ■? . Mb , o o , 



where r = Mj^(*) / Mb- Putting the light (anti-)quark momenta in P, , mesons as 
/ci, A;2, and /cs, respectively, we can choose 



ki = (xiP+, 0, kir), k2 = (X2P2+, 0, k2T), ks = (0, X3P3-, kgr). 



(4) 
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Unlike QCD factorization approach, we do not neglect the transverse momentum kj^ in 
the above expressions, by which to avoid the endpoint singularity. 

If the decay involves one or two vector mesons in final states, the longitudinal polar- 
ization vectors eg , up to the order of are given by 



and transverse polarization vectors ef", ej are 

e^=(0,0,lT), e[ = (0,0,lT). (6) 
Then, integration over Jv-^ , 7 Biiici kn in eq. (0) leads to: 



Amplitude ~ J dxidx2dxzhidhib2dh2bzdhz 

Tr 6i)$d(x2, &2)$d(x3, h)H{x,, k, t)St{xi) 6"^^] , (7) 

where bi is the conjugate space coordinate of the transverse momentum kiT, and t is the 
largest energy scale in H, as a function in terms of Xj and fej. The last term, e"'^*^*-*, contains 
two kinds of logarithms. One of the large logarithms is due to the renormalization of ultra- 
violet divergence In tb, the other is double logarithm In^ b from the overlap of coUinear and 
soft gluon corrections. This Sudakov form factor suppresses the soft dynamics effectively 
ll| , so it makes a perturbative calculation of hard part H applicable at the intermediate 
scale. 

As a heavy meson, the B meson wave function is not well defined, so is D meson. In 
heavy quark limit, we may use only one independent distribution amplitude for each of 
them. 

6) = [/> + Mb] i,<Pb{x, 6), (8) 

b) = -^75 [/' + Md] 0d(x, b). (9) 
For the vector D* meson, it is expressed as: 

^D' (X, 6) = + Mz).] (x, b) . (10) 

The hard part H, which is channel dependent, can be calculated perturbatively. We 
show the calculated formulas below for different channels. 
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A. D+D-, 5° ^ D+D- decays 

In the decay B'^ DfD~, the effective Hamiltonian at scale lower than Mw is jl^ : 

= ^'^:iycd [CMO^i^i) + CM02{^^)] , (11) 

Oi = {db)v-A{cc)v-A, O2 = {cb)v-A{dc)v-A- (12) 

In above functions, Ci^2{iJ) are Wilson coefficients at renormalization scale /x. And sum- 
ming over SU(3)c color's index a, QaYi'^—lblQai are abbreviated to {qq')v-A- Penguin 
operators may also have contribution, but they usually have smaller Wilson coefficients. 
Here we neglect these diagrams. The lowest order diagrams for the hard part H calcu- 
lation, are drawn in FigQ according to this effective Hamiltonian. Just as what we said 
above, there are only annihilation diagrams. 

For the decay -8° D^D~ , the effective Hamiltonian at scale lower than Mw is |l^ : 

^eff = [c,{^^)0[{^^) + C2{M{^^)] , (is) 

0[ = isb)v^Aicc)v^A, = icb)v-Aisc)v-A- (14) 

Comparing with Eqs. pi|12|l . the only changes in Eqs. pH|14|l are the replacements of the 
CKM factor Vd —>■ Vcs and the quark d s. As we will see later the branching ratio of 
D^D^ will be much larger than that of B^ D+D^ decay, because of this larger 
CKM factor Vcs- The lowest order diagrams for the hard part H calculation, are then 
similar to B^ decay in Fig^only replacing the d quark by s quark. 

In decay B^ DfD~, we get the following analytic formulas by calculating the hard 
part H at first order in a^. The factorizable annihilation diagrams in FigCJi and b cancels 
each other, which is a result of conservation of vector current and parity invariance. 

With the meson wave functions, the decay amplitude for the nonfactorizable annihila- 
tion diagrams in FiglH^c) and (d) results in 
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Ma = —^GAttCfM'^ / dxidx2dx3 i 
v6 Jq Jo 



bidbi b2db2 <Pb{xi, bi)4>D,{x2, &2)0D,(a;3, 



-X'iE^{t]^)h^^\xi,X2,X'iM,h) + X2E^{tl^)hf\xi,X2,X2„bi,b2) , (15) 

where Cp = 4/3 is the group factor of SU(3)c gauge group. The function Em is defined 
as 

Era{t) = C2{t)as{t) e-Ss{t)-5z,(t)-5z,W^ (^g) 

where Sb, Sd result from summing double logarithms caused by infrared gluon corrections 
and single logarithms due to the renormalization of ultra-violet divergence j^. 
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The functions h^a^ and h^a^ are the Fourier transformation of virtual quark and gluon 
propagators. They are defined by 

where Hq^^(z) = Jo(-2) + iYo{z), and -F(j)S are defined by 

(1) 

F{2) =xi+ X2XS + (1 - r'^){x2 + xs- XiXs - 2x2X-i). (19) 

The hard scale t's in the amplitudes are taken as the largest energy scale in the H to kill 
the large logarithmic radiative corrections: 



F?^) = XiXsi^l - r^) - X2X3{1 - 2r^), 



X2X3A/biA/h)- (20) 



(i) 

Applying the power counting rule established in ref.0, 0], we keep only the leading 
order contribution of r expansion in the numerator of the above equation p7|) . The 
hierarchy relation Aqcd itt-dm <^ is assumed. The terms are kept in the 
denominators of (|17p. since it may sometimes affect the imaginary part heavily. It is 
easy to see that the momentum carried by the intermediate gluon is Mb \/x2Xs{l — 2r^) 
which is only suppressed by a factor of a/(1 — 2r^), comparing with that of the charmless 
B ^ mr decay :2j. In the heavy quark limit, r — 0, the momentum of the gluon is the 
same for the two kinds of decays. The formulas derived here support the argument at the 
introduction that perturbative calculation is still applicable to the B DD decays. 

The decay width F for B^ —>■ DfD~ decay is then given by 

r{B' DtD-) = ^^(1 - 2r')\v:,V^,M^\'. (21) 

i/OTT 

Similar to decay B^ D^D^ , the width for B^ D^D^ is 

F(i?° ^ D^D-) = ^^(1 - 2r2)|\/;KA|^ (22) 

One need only replace the Dg wave function (j)£,^ by D meson. Enhanced by the CKM 
factor |V^|,KsP, the B^ decay width will be larger than that of 5° decay. 
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B. B° D*^D-, D+D*- and 5° ^ D*+ , D+D*- decays 

For final states with one pseudo-scalar and one vector mesons, only the longitudinal 
polarization of vector meson contribute. The decay amplitude takes the same form as 
the amplitude of B to two pseudo-scalar mesons (fT3j) . For decay —>■ D*~^Dj, one 
need only replace one of the Dg meson distribution amplitude by D* one. The width of 
B^ DfD*~ must have the same width as B^ D*'^Dj . The contributions of Figl.(a) 
and (b) can not be cancelled by each other because of the difference of fc^ and fn*, but 
it is still negligible. 

Accordingly, the B^ — > D*^D^ , D^D*^ decay amplitudes also take the same form as 
B2 D+D-. 



C. 5° ^ D*+D*- and 5° ^ D*+D*- decays 

There are contributions not only from the longitudinal polarization but also from two 
transverse polarizations in i? ^ VV decays, where V denotes the vector meson. Therefore 
the decays 5^ D*^D*^ and B^ D*^D*^ are more complicated than B PP or 
B PV. 

In the covariant form, the decay amplitudes of non-factorizable annihilation diagrams 

are 

pi poo 

M'^ = ^=128TTr'^CF / dxidx2dx3 / bidbib2db2 4>B{xi,bi)(f)D*{^2,b2)4>Dt{x3,b2) 
y^i^c Jo Jo 

X |(e2 ■ e3){p3- P2){X2+ Xs) E^{t]^)h^^\xi,X2,X2.,bi,b2) - Em{tUh'^^\xi,X2,X3,bi,b2 

+ (e2 ■P3)(e3 ■P2)(a;2 -3:3) E^{tl^)h''^\xi,X2,X3,bi,b2) + E„(t^)/ii^^(a;i, 0:2, 0:3, 61, 62) 

-ie'"'^'^e2f,e3uP2pP3aix2-X3) Em(tljh^^\xi,X2, X3, bi, b2)-Em{tln)h'^a\xi, X2, X3, bi, 62 



(23) 

with the convention tr(75 P ^ fl) = —Aie^^^'^aabpc-ydp and e°^^^ = 1. Just as Section 
pi Aj) . the contributions of diagrams(a) and (b) cancel each other. 

If we set the to be longitudinal polarization only, the above formula goes back 
to the eq. ljlSp . From the above functions, we can also see that the contributions of 
transverse polarizations are proportional to factors of r^, which are suppressed comparing 

with longitudinal ones. In our calculation, we set rric ~ just because m^(*) — mc ~ 

(3) (3) 

A. And A/m (,) — in the heavy quark limit. 
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III. NUMERICAL RESULTS 



For B meson, we use the same wave functions as other charmless B decays 0, 0| , which 
is chosen as 



')b{x, h) = Nbx^{1 — xY exp 



Ml 



(24) 



The parameters = 0.4 GeV, and Nb = 91.745GeV which is the normahzation constant 
using /b = 190MeV, are constrained by charmless B decays J HI. For Bg meson, we 
use the same wave function according to SU(3) symmetry. That is Ub = 0.4 GeV, but 
Nb, = 119.4GeV, using /b, = 236MeV. 

For D(^gy the distribution amplitude is taken as 



/ {.)a;(l - a;){l + a 



D 



2x)}. 



(25) 



Since the heavy wave function is less constrained, we use a^ic*) = 0.6 ~ 0.8 GeV and 
a„{,) = 0.2 ~ 0.4 GeV to explore the sensitivity of parameters. Other parameters, such 

^ s 

as meson mass, decay constants, the CKM matrix elements and the lifetime of 5(5) meson 
[lit are given in Table HI 



TABLE L Parameters we used in numerical calculation. 





m^o = 5.28GeV 




= 5.37GeV 


Mass 


m£,± = 1.87GeV 




= 2.01GeV 




mj^± = 1.97GeV 




= 2.11GeV 


Decay 


= 240 MeV 


fD*± 


= 230MeV 


Constants 


= 241MeV 




= 211 MeV 


CKM 


l^cdl =0.224 


\Vcb 


= 0.042 




\Vcs\ = 0.974 






Lifetime 


tbo = L54 X 10-12 s 


tbo = 1.46 X IQ-^'^s 



The calculated branching ratios in PQCD are sensitive to various parameters, such as 
the parameters in wave functions of B and D^ly Because the wave functions are from 
non-perturbative effect, we can not define them exactly. In Table IIT| we show examples 
of the sensitivity of the branching ratios to parameters in ()25|1 . The predictions of PQCD 
depend heavily on and a^i, which characterize the shape of D) ! wave function. 
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TABLE II: The sensitivity of the decay branching ratios to change of a_D(aDs)- 



a-Ds 


Br(SO - 




-)(io-5) 


0.2 




6.2 




0.3 




7.8 




0.4 




9.8 




an 


Br(5° - 


D+D- 


-)(io-3) 


0.6 




3.3 




0.7 




3.9 




0.8 




4.5 





From above discussion, the branching ratios within the reasonable range of parameters 
in wave functions are given as 

Br{B^ DtD-) = (7.8±?:|]) x 10-^ Br(50 ^ D+D-) = (3.6 ± 0.6) x IQ-^; 
Br(50 D*+D-) = {6.0±\i) x 10-^ Br(5° D*+D-) = (3.6 ± 0.6) x IQ-^; 
Br(5° D*+D*-) = {8.5±li) x 10~^ Br(5° ^ D*+D*-) = (6.4 ± 1.2) x lO'^. 

(26) 

Because Jd* is smaller than f^^, we can see the branch ratio of 5° — D*~^D~ is a little 
smaller than D^D~. The Br(i?° — > D*~^D*~) is larger than the others because of 

the extra contribution from transverse polarization. 

In the decay 5 — » vrvr, m.^ is much lighter than m^. The energy release in the decay 
is very large. The final tt mesons runs very fast and they may not have enough time to 
exchange the soft gluons and resonance. Recently, the B — > Dn, B — > LKK{(f)) decays 
with one heavy Z?|*j meson in final state are calculated in PQCD approach 0,E|- And the 
results are consistent with the experiments, which shows that the final state interaction 
may not be important in those decays, although the energy release is smaller than that 
in 5 ^ TTvr decays. Here we also calculate B^ and 5° £)(*)+£)(*)- 

decays with two heavy mesons in final states in PQCD. We get large branching ratios 
comparable to other predictions j^. This may be a hint for PQCD to work good for these 
decays. The soft final state interaction in those decays, for example, B^ D^D~, and 
D^D^ D^D^ through exchanging K^{K^) is somehow smaller than the perturbative 
picture. 

For consistent check, in Figure |2l we show the contribution to the branching ratio 
of decay B^ — > DfD~ from different ranges of Os/tt, where the hard scale t is given 
in appendix. From this figure, we find that most of contribution comes from the range 
Os/tt < 0.3, implying that the average scale is around Aqcd^b- It is then numerically 
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FIG. 2: Contributions to the branching ratio of decay B — > D^D^ from different ranges of 
as/n. 



confirmed that PQCD may be even applicable to 5° DT'^Dr and 5° !:)(*)+£)(*)- 
decays. 

In ref. P, J.O. Eeg et al. computed those decays using heavy-light Chiral quark model 
and their results read: 

Bi{B^ DfD-) = 7.0 X 10-^ Br(50 ^ D+D-) = 1.0 x lO'^. (27) 

Obviously, for decay 5° D^D^, we have the same result. But we got different result 
for decay D^D^ though our results are at the same order. Unfortunately, there 

is no direct experimental result about these decays up to now. We hope those branching 
ratios will be measured soon in future and these two theories can be tested. 



IV. CONCLUSION 

In this paper, we try to estimate the branching ratios of Di*^'^Di*^ and 

£){*)+ jji*)- decays in the heavy quark limit using perturbative QCD approach. These 
decays can occur only through annihilation diagrams because the four quarks in the final 
states are not the same as the ones in B meson. Our numerical results agree with the 
heavy-light chiral quark model for B^ D+D^ and 5° D^D^ decays, which is not 
very small. We also give large branching ratios for channels with one or two vector mesons 
in final states. There is a hint for not large soft final state interactions. We hope new 
experimental results will give a test for our results. 
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